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t I . A complete and rigorous determination of the possible ground states for D-wave pairing Bose 

' condensates is presented. Using an orbit space approach to the problem, we find 15 allowed phases 

(besides the unbroken one), with different symmetries, that we thoroughly determine, specifying the 
group-subgroup relations between bordering phases. 
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I. INTRODUCTION 



^ i' Superfluidity and superconductivity are justified on the basis of the macroscopic condensation of Bose quasi- 
^ ' particles. The classical Bardeen, Cooper and SchriefFer theory for superconductivity dates 1957. Soon after, a BCS- 
type transition was proposed for the Fermi system '^He by Anderson and Morel |l| . Cooper pair formation was thought 
' to occur in an L state, to take into account the hard core nature of '^He atoms interaction. The superfluid phases 
were actually observed and the nature of p-wave pairing is now well established for '^He. The theory oi L ^ 0- 
superfluids is also relevant for "unconventional" superconductivity. The high temperature superconducting (HTS) 
'■^ oxides are anomalous in their non-Fermi liquid normal state properties and share with heavy fermion superconductors 
^ unconventional d-wave pairing 1^,^. The underlying microscopic mechanism inducing superconductivity in these 
O ' materials is still unclear and is one of today's major challenges [^,||. 

I ^ I ■ Such a situation has motivated the efforts at studying the macroscopic properties of unconventional superconductors 
through the Landau theory of phase transitions [Q. Moreover, the identification of the order parameter symmetry 
may be considered as a preliminary task in the construction of viable models describing the attractive nature of the 
^ pairing interaction. 

Direct experimental evidence of an order parameter unconventional structure lies on multiple phase transitions. The 
heavy fermion compounds Ui-xThjjBeis display four different superconducting phases in the T-x phase diagram and 
[ UPta displays three superconducting phases in a T-H phase diagram 1^. Furthermore, from power-law temperature 
' behaviour of thermodynamic and transport properties (e.g. specific heat, magnetic penetration depth), a non-trivial 
' node structure for the gap function may be inferred, compatible with high i-pairing. That is the case also for HTS 
t oxides, for which, at present, clear proofs are lacking of the existence of more than one superconducting state. 

Actually, it must be pointed out that, unlike superfluid He, which is an isotropic fermion system, Bloch electrons 
in a superconductor crystal lattice exhibit, in general, a reduced finite symmetry; in addition, the existence of 
y imperfections, even in the cleanest samples, can partly destroy the gap node structure at low energy. So, the very 
fact that the experimental observation cannot at present completely unravel the node structure of the gap function, 
^ reinforces the necessity of classifying all the possible symmetry breaking schemes. 

Q _ Moreover, the role of Fermi isotropic space beyond a zero-th order approximation in phenomenological theories, is 
O . still supported by some recent studies: In HTS oxides, low symmetrical crystal fields (tetragonal and orthorombic), 
^ ' as well as spontaneous strain, have weak infiuence on the temperature of the superconducting phase transition ||^ or 

• i-H , on the penetration depth ||^ . This also means that low symmetry crystal fields do not directly infiuence condensate 

X '■ states. 

^ ' The possible ground states of a high L superfluid has been the object ofjntense investigations during the 60's and 



ICj] identified five different 
rtll analyzed the problem 



the 70's. ^From the solution of the state equations, Anderson and Morel ||^ and Mermin 
phases, through the minimization of a 4-th degree Landau potential. Schakel and Bias 

using only group theoretical arguments. Capel and Schakel |l^, taking advantage of the results in jlll, computed 
the possible ground states of condensates driven by p- waves. They also investigated the consequences of lowering 
the residual symmetry resulting from the solutions of this problem by means of strong spin-orbit interactions. Since, 
in this case, the Cooper pair is in a J = 2 state {S = 1,L = 1), the order parameter is represented by a traceless 
symmetric tensor. The same holds true for an {S ~ 0, L = 2) state order parameter, so those results may be directly 
applied to D-wave pairing. According to the authors of ref. [O], eleven different phases are allowed. Their analysis. 
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however, is questionable, since time reversal symmetry cannot be neglected in the theory of condensate states and 
the use of a 4-th degree polynomial Landau free energy strongly limits the number of phases that would be allowed 
by the symmetry of the system. 

Our aim in this paper is to give a definitive answer to the classification of possible symmetry breaking patterns in 
D-wave pairing Bose condensate, in the framework of the Landau theory of phase transitions. To this end we make 
use of a geometrical invariant theory approach to the problem, proposed in ref. The strategy is to exploit a set 
of basic invariant polynomials of the symmetry group of the system, as fundamental variables in the description of 
the phase space of the system and in the minimization procedure of the free energy . 



II. THE ORBIT SPACE APPROACH 

Let us briefly recall some basic elements of the orbit space approach To this end, we shall denote by x € R" 
a vector order parameter, transforming linearly and orthogonally under the compact real symmetry group G, and by 
$(a; x) the G-invariant free energy, expressed also in terms of state variables a. The points xo(a), where the function 
4>a{x) — $(a;a;) takes on its absolute minimum, correspond to the stable phase of the system, whose symmetry is 
determined by the isotropy subgroup, Gxa , oiG aX xq. Owing to G-invariance, the stationary points of the free energy 
are degenerate along G-orbits. Since the isotropy subgroups of G at points of the same orbit are conjugate in G, 
only the conjugate class, [G^o], of Gxa ^-S- the symmetry (or orbit-type) of the orbit through xq, is physically 

relevant. 

The set of all G-orbits, endowed with the quotient topology and differentiable structure, forms the orbit space, 
R"/G, of G and the subset of all the G-orbits with a given symmetry forms a stratum of R"/G. Phase transitions 
take place when, by varying the values of the a's, the point xo(a) is shifted to an orbit lying on a different stratum. 
If $(a;a;) is a sufficiently general function of the a's, by varying these parameters, it is possible to shift a;o(a) on 
whichever stratum of R"/G. So, the strata are in one-to-one correspondence with the symmetry phases allowed by 
the G-invariance of the free energy. On the contrary, extra restrictions on the form of the free energy function, not 
coming from G-symmetry requirements (e.g., the assumption that the free energy is a polynomial of low degree), can 
limit the number of allowed phases. 

Being constant along each G-orbit, the free energy may be conveniently thought of as a function defined in the 
orbit space of G. This fact can be formalized using some basic results of invariant theory. In fact, the G-invariant 
polynomial functions separate the G-orbits, meaning that, for any two distinct G-orbits, there is at least a polynomial 
G-invariant function assuming different values on them. Moreover, every G-invariant polynomials can be built as real 
polynomial functions of a finite set, {pi{x), . . . ,Pq(a;)}, of basic polynomial invariants [integrity basis of the ring of 
G-invariant polynomials), which need not, for general compact groups, be algebraically independent. The number 
of algebraically independent elements in a em minimal set of basic polynomial invariants is n — v, where v is the 
dimension of the generic (principal) orbits of G. Information on the number and degrees of a minimal set of basic 
invariants, and the degrees of the algebraic relations [syzygies) among them, can be inferred from the Molien function 
of G. 

Let us call go the number of independent elements of the set {p}. The range of the orbit map, x i-^ p{x) = 
{pi{x), . . . ,pq{x)) e R"?, yields a realization of the orbit space of the linear group G, as a connected semi-algebraic 
surface, i.e. a subset of R', determined by algebraic equations and inequalities. The orbit space of G is, therefore, 
a closed and connected region of a go-dimensional algebraic surface, delimited by lower dimensional semi-algebraic 
surfaces. 

Like all semi-algebraic sets, the orbit space of G presents a natural .stratification, since it can be considered as the 
disjoint union of semi-algebraic subsets of various dimensions [geometrical strata), each stratum being in the border of 
a higher dimensional stratum, but for the highest dimensional one, which is unique [principal stratum). The connected 
components of the symmetry strata are in one-to-one correspondence with the geometrical strata. The symmetries 
of two bordering strata are related by a group subgroup relation and the lower dimensional stratum has a larger 
symmetry. 

The orbit space can be identified with the semi- algebraic variety, S, formed by the points p G R*, satisfying the 
following conditions i) and ii) [^3|jl^]: 

i) p lies on the surface, Z, defined by the syzygies; 

ii) the q X q matrix P[p), defined by the relation 



n 

Pab[p[x)) = ^ djPa[x) djPb[x), Va; e R' 



(1) 



is positive seinidefinite at p. 

The relations defining the strata can be obtained as positivity and rank conditions on the matrix P{p) and the 
minimum of '^{a]x) can be computed as a constrained minimum in orbit space of the function $(a;p), 

$(a;p(a;)) = $(a;x), Va; € R" (2) 

or from the solutions of the equation jis) 

q 

Y,Pab{p)db^{a-p)^Q, a=l,...,g, (3) 

b=l 

which is equivalent to the state equation, 9$(a; x) / dxj = 0, j — 1, . . . ,n. 



III. SYMMETRY OF THE ALLOWED D-WAVE CONDENSATE STATES IN ISOTROPIC SPACE 

The formation of D-wave condensate states breaks the symmetry of the isotropic 3-dimensional space, which corre- 
sponds to the group 03(8)Ui x (T), where O3 is the complete rotation group, Ui is the group of gauge transformations 
and (T) is the group generated by the time reversal operator T. 

The symmetry of the allowed D-wave condensate ground states is defined by the relative values of the complex 
coefficients in the decomposition of the gap-function, A, in terms of spherical harmonics with L = 2: 

2 

A(0,0)= D„^Yr{e,<P) (4) 

m=-2 

The set of functions {Y™ ,Y™*} yields a basis of a ten-dimensional (10 D) space hosting a real representation of 
the symmetry group 03(8)Ui x (T). A general element, 7, of the group will be denoted by a triple 7 = {p, e'^ , e), 
where, p S O3, < </> < 27r and e ~ —1, or -1-1 according as a time reflection is involved in the transformation, or not. 

The action of G can be transferred to a real irreducible action on the 10 D vector formed by the coefficients 
{D2, . . . , D_2^ • • • 1 ^2}- The representation of G thus obtained can be realized in the 10 D real vector space of 

a couple of two independent, real, second rank, symmetric, traceless tensors, x[\j'^ and x\^\ i,j — 1,2,3, which can 
be considered as the real and imaginary parts of a complex 3x3 matrix ?/;, whose elements will be written in terms 
of five complex coordinates, Zj: 

Zj=Xj+ix5+j, j = l,...,5; XiGR, (5) 
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Z3 



-Zl 



15 

73 



Z2 




The coordinates Da are connected to the Zj by the following relations: 



iZi+ Z2 



Di = 



i Z3 -I- Z4 
^2~' 



Do = Z5, D-i = 



IZ3- Zi 



V2 



D-2 = 



IZi — Z2 



V2 



The matrix ip transforms in the following way under a general transformation 7 = (p, 0, e) S G: 



7 • ip 



p 



(6) 



(7) 



(8) 



where ip' = ^ or according as e = -|-1 or —1 and the apex T denotes transposition. As a consequence, the group 
G acts as a group of linear, real, orthogonal transformations on the vector order parameter x G R^". 

The kernel of the representation of G just defined is the group generated by the space reflection. So, it will not be 
restrictive to assume that the symmetry group is G = SO3 (g) Ui x (T) and, when speaking of G, in the following, we 
shall always refer to this linear group acting in the vector space R^°. 
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The linear group G has a trivial principal isotropy subgroup (the isotropy subgroup of generic points of R"*^*^), 
thus the principal G-orbits have the same dimensions as G and its orbit space, i.e., the quotient space 'R}^/G, has 
dimensions go = 10 — 4 = 6. 

The Molien function of G, M{rj), can be calculated in the form of an invariant Haar integral over G (see, for 
instance, 

Mi,) = I 1^1 < 1' (9) 

Jg det(l-?7g) 

where ^{g) is a normalized invariant measure on the group G, the integration is over the whole group G and g ^ G. 
An explicit calculation of the integral leads to 

(l-?72)(l_j^4)2(l_^6)2(i_^8)- ^ > 

Equation ( |lo|) yields the following indications, whose validity has been checked through direct calculations: 

1. A minimal integrity basis for the linear group G contains nine elements, {pi, . . . ,pq\ with degrees {di, . . . , dg) — 
(2,4,4,6,6,8,8,10,12). 

2. The invariants pi are connected by five independent syzygies of degrees 16, 18, 20, 22 and 24. 

3. The most general G-invariant polynomial, like a general non-equilibrium polynomial Landau potential, ^{a;p), 
can be written as a polynomial function of the elements of the integrity basis {pi}i=i...9, in terms of five arbitrary 
polynomials Qi = Qi{a;pi, . . . ,pq), j = 0, . . . ,4: 



The elements of the minimal integrity basis can be chosen in the following form: 

Tr(V'2)Tr(V'2V*)Tr(V'*^)] , 
Tr(V'*2)(Tr(V'V*))'] , 
Tr(V'2V'*)Tr(V'^) (Tr(V'*2))' 



= Tr(V^V*)-n=i2;| 

P2 = Tr[(^V*)2] , 

P3 = |Tr(V'2)|2, 

P4 = |Tr(V3)|2, 

p, = |Tr(V'V)l', 



p& 


= ^ 


P7 


= 3? 


P8 




P9 


= 3? 



(12) 



(Tr(v^2))' {TTir')y 



Using these definitions, we have determined the explicit form of the syzygies, of the P-matrix elements, of the 
equations and inequalities determining the strata in the orbit space of G. For each stratum, denoted by S'--'^''^\ where 
d denotes the dimension and r is an enumeration index, we have picked up a "typical point" and determined the 
corresponding isotropy subgroup of G. All these results are essential pre-requisites for a rigorous calculation of the 
minima of a G-invariant polynomial along the lines indicated in the Introduction. 

Part of our results are resumed in Tables |, || and III and in Fig. 1, where the group-subgroup relations among 
the symmetries of bordering strata are also specified. 

The explicit expressions of the syzygies and of the elements of the P{p) matrix and the relations defining higher 
dimensional strata would require too much space to be written down here. 

As a simple example of the effectiveness of our approach, we have also repeated Mermin's calculation of the minimum 
of a general fourth degree polynomial G-invariant free energy 



$W(p)=a„|-f^a, 



Pv 



(13) 



in the additional assumptions that it is bounded below and has a local maximum at the origin (ai < 0). The cti's are 
connected to the parameters a, Pi, P2, Pa used by Mermin by the following relations: 



ao = 2/?2 + /33 



ai = a, 



a2 



2 



"3 = /9i + 



(14) 



With the definitions: 
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P^'"~d72' P = (P2, ■ • • ,P9), (15) 

Pi 

and 

A(p) = q;o + 2a2P2 + 2a3P3, (16) 
the polynomial (p) can be put in the following convenient form: 

$(4)(p)^^A(p) + aiPi. (17) 

Since, owing to its definition, pi ranges over the whole non negative real numbers, ^^^\p)\p^g is bounded below 
(stability condition) if and only if the minimum, S, of A(p)|p^_j^ is positive. Being the minimum of the r.h.s. of ([T^), 
thought of as a function only of pi > 0, equal to — q;i/(2A), the absolute minimum of (p) [^^^ is —af/ (25). In this 

way, the problems of rendering explicit the stability condition and evaluating the minimum of *J'(?')|pg5 ^^'6 reduced 
to the calculation of 6. 

The absolute minimum of 6 in each singular stratum can be easily obtained from the equations of the strata, listed 
in Tables || and |l| for strata with dimensions < 3. For the principal stratum, it is easier to solve the projection of 
equation (^ in the unit sphere (X]i=i ^^f = Pi = 1)^ which, using (15), can be written in the form | |2l| 



J=2 



dAjp ) 
dpj 



y(p)^F^-0 (18) 



and to select, subsequently, the solutions lying in the principal stratum. 

A comparison of the values of the minima in the different strata, obtained in this way, leads to the results resumed 
in Table IV and illustrated in Figure 2. Owing to the low degree of the polynomial defining $*^''-'(p) in (|l3|), and 



the consequent low number of free parameters a, the absolute minimum presents strong degeneracy, particularly for 
special values of the a's. If these special values are excluded, spontaneous breaking of the symmetry can generate 
only five distinct phases out of fifteen permitted by the G-symmetry; some of them are unstable. For no non trivial 
values of (a2, as) does the absolute minimum he on the stratum 5^^'^'. 

For general values of the a's, our results are in agreement with Mermin's ones. Let us add a few words about 



the perturbative stability of the three degenerate phases in the region TZ2 (see Table IV). For (00,02,0:3) S TZ2, 
the addition to the free energy, $'^^^(p), of a "small" perturbation, consisting in an invariant polynomial of degree 

SIX, 

e(6) 

= 04P4 + 05P5, sphts the three degenerate minima determined by the 4-th degreee termQ This is easy to 
check, at least in the additional assumption that the perturbation leaves the absolute minimum in one of the strata 
corresponding to the degenerate phases. In fact, at the first perturbative order, one obtains from Tables || and|l| the 
following shifts, Q^^^y in the values of the 6-th order free energy at the points where ^^^\p) takes on its degenerate 
absolute minimum under consideration: 

QS) = (t)'^' QS4) = (^)^"4 + «a)C, (19) 

where < ^ < 1/6. 

Since —ai/6 > 0, the absolute minimum will be perturbatively stable on 5'^'^^ or, respectively, S'--^'^\ according as 
0(4 + as is negative or positive. 

As stressed in the Introduction, the difficulties mentioned above can be overcome if one puts less restrictive upper 
limits to the degree of the polynomial describing the free energy. It is trivial, for instance, to realize that the following 
class of bounded below polynomial functions have a vanishing maximum at the origin of R^" and display an 
absolute minimum at the arbitrarily chosen point p (z S: 
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[(P.-P.)'"'-P^], (20) 



^The inclusion of terms of 6-th degree, depending only on pi,P2 and ps would be useless for splitting the degenerate minima, 
so these terms will be neglected, with no loss of generality. 
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where the a's arc positive constants and the n's are positive integers. 

The physical imphcations of our results and the derivation of a more realistic form of the free energy will be discussed 
in forthcoming papers. 
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TABLE L Relations defining strata, S^^'''', of dimensions 1 in orbit space. For 2 < i < 9, g, = Pi/(pf*^^) and dj denotes the 
degree of the polynomial Pi{x) (x is the order parameter). 
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1/6 
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TABLE II. Relations defining strata, S^^'^\ of dimensions 2 in orbit space. The pi's are defined as in Table | and e = ±1. 



p\r 1 2 3 4 5 

~T2 (2+e^)/6 172 {2+e)/G 172 r~ 

P3 e 12-2^ 

P4 (2-0^(1 + 0/12 (2-0^(1+0/12 ^ 

P6 0(1+0/12 ^ 

P6 (2 -o(i + ocVi2 e 

P7 (i + o^Vi2 e 

PS (2-0(l + 0?Vl2 e 

P9 (2-0'(l + OCVl2 ^ 

grange ]0,i^[ ]0,1[ ]0,e[ ]0, i[ ]i,l[ 



TABLE III. Possible symmetry strata, S'*'*'''^ for D -wave driven pairing in isotropic space (d denotes the dimension of 
the stratum and r is a secondary enumeration index). From left to right, the columns refer to phase reference numbers 
according to our((ci, r)) and, respectively, ref.[12] (A/") classification, order {\H\) and residual symmetry group (H), complex 
coordinates {zj — Xj + ix^+j, j = 1, ... ,5) of the //-invariant superconducting vector order parameter and corresponding 
values of the partial wave amplitudes (D2,...,d_2)- In columns five and six, the t's are real coordinates, while the v's are complex 
ones. The notations are the same as in [19], Rz{(t>) denotes a clockwise proper rotation of the vectors, about the z-axis and 

O2 = {Rz{(f')}o<4,<2-^ U {C2i:R2(';/))}o<<;,<27r- 



(d,r) 


JV 


\H\ 


H 


(21,..., Z5) 


V2{D2,Di,Do,D-i,D^2) 


(1,1) 


II 


OQ 


(C2,r>x{R,(«Ui(2</,)}^ 


(it, -t, 0,0,0) 


{2t, 0,0, 0,0) 


(1,2) 


VIII 


00 




(0,0,0,0,4) 


(0,0,tV2,0,0) 


(1,3) 


IX,X 


16 


(C2^,r,C4,Ui(7r)) 


(0,t, 0,0,0) 


{-t,0,0,0,-t) 


(1,4) 


XI 


24 


(C2x,C2aT,C3jUi(47r/3)) 


(0, -it, 0,0,4) 


(it,0,tv'2,0,it) 


(1,5) 


I 


00 




(0,0,t,-it,0) 


(0,0,0,2 it,0) 


(2,1) 


V 


6 


(C2x-r,C3,Ui(47r/3)) 


(iti, — ti,t2,— it2,0) 


(ti,0,0,it2,0) 


(2,2) 


IV 


4 


<C2xT,C2,Ui(7r)) 


(0,0,ti,it2,0) 


(0,iti+i t2,0,iti— it2) 


(2,3) 
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(0,iti,0,0,t2) 


(-iti,0,t2\/2,0,-iti) 


(2,4) 


VII 


8 


{C2^,C2,,r) 


(0,ti,0,0,t2) 


(-ti,0,t2y2,0,-ti) 


(2,5) 


VI 


8 


(C2xr,C4,Ui(7r)) 


(iti,t2, 0,0,0) 


(ti~t2,0,0,0,-ti-t2) 


(3,1) 


III 


4 


(C2^,C2xT> 


(iti,t2,0,0,t3) 


(ti-t2,0,t3 V2,0,-ti-t2) 


(3,2) 




4 


(C^2x,C^2s) 


(0,ui ,0,0,1)3) 


(-1)1,0, 1)2 \^,0,-Di) 


(4,1) 




2 


<C2xT> 


(iti,t2,t3,it4,t5) 


(ti — 12 ,it3+it4,t5y2,it3— it4, — ti — 12) 


(4,2) 




2 


{C2.) 


(iJl, -1)2 ,0,0, 1)3) 


( — iDl— 1)2, 0,^3 V^,0,iDi — 1)2 ) 


(6,1) 




1 


{1} 


(ui, 1,2,113,1)4,1)5) 


( — iDl — D2,iD3-t-i'4,D5 N/2,iD3 — V4 ,ivi —V2) 



TABLE IV. Absolute minimum, $ 



(4) 



— Qfi/(25), of a general, bounded below, G- invariant 4-th degree polynomial, 
"l?'-'*^ (Qf,p) = QfoPi/2 + X/j=i'^jPj' "^1 0, and hosting strata, S'-'^'^\ as functions of the coefficients a. The denomination 
of the strata is the same as in Table |ll|. 



a range 



(d,r) 



7^l : Max(0, -3q:o/2, -603) < Q2 
TZ2 : —603 > 02 > Max(— ao — 2a3, 203) 
7^3 : -0:0/2 < Q2 < Min(0,2Q:3) 
7?.i3 : = 02 < Min(ao,a3) 
7^l2 ; Max(-3ao/2,0) < 02 = -603 
72.23 : —ao/2 < a2 — 2as < 
7?.i23 : a2 = Q3 = < ao 



ao + 2a2/3 
ao + a2 + 2a3 
ao + 2a2 

ao 
ao — 4a3 
ao + 4a3 
ao 



(1,4) 
(1,2), (1,3), (2,4) 

(1,1) 

(1.1) , (1,4), (1,5), (2,1), (3,1), (4,1) 

(1.2) , (1,3), (1,4), (2,3), (2,4), (3,2) 
(1,1), (1,2), (1,3), (2,4), (2,5) 

all, except (0,1) 
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FIG. 1. Possible phase transitions between bordering strata, connected, in the figure, by continuous sequences of one or more 
arrows. The notations are the same as in Table III. 



S<°) = {0} 
SO3 (g) Ui X T 




5(4,1) 




S(4,2) 









FIG. 2. Localization of the absolute minimum of a fourth degree G- invariant polynomial, (a,p) ~ aoPi/2 + ^ ctjPj, 
Qfi < 0, as a function of its coefficients. For values of (ao,«2,a3) in TZ\ or in TZ^ or in TZs, the absolute minimum lies, 
respectively, in the strata S^^'^'^ or {S^^'^' , S^^'^^ S'^'"''} (degenerate minimum) or S'^'^-*. For particular values of the a's, see 
Table livl 



8 




9 



